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FOR ORBIT INSIDE EARTH'S 
MAGNETOSPHERE, START 
HERE: 



FOR ORBITS AT OR BEYOND 
GEOSYNCHRONOUS, START 
HERE: 



AUXILIARY PROGRAMS YOU 
MAY WISH TO RUN: 



GTRN: 

EVALUATES 
GEOMAGNETIC SHIELDING 



200 



FLUX: 

EVALUATES FLUXES AT EXTERNAL 
SURFACE OF SPACECRAFT 



300- 



■100 



TRANS: 

TRANSPORTS NUCLEI 
THROUGH SHIELDING 



PUP: 

PROTON-INDUCED 
SEU RATES 



T 

500 




SHIELDFILE: 

FORMATES USER-SUPPLIED 
SHIELDING DISTRIBUTION FOR 
INPUT TO TRANS 



400 



CALCULATES 
LINEAR ENERGY 
TRANSFER (LET) 
SPECTRA 



•600 



HUP: 

DIRECT-IONIZATION 
INDUCED SEU 
RATES 



-700 



FIG. 1 
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INPUT PARTICLE RIGIDITY IN GV, 
ARRIVAL DIRECTIONS (I.E. ZENITH 
ANGLES e m in&e max AND azimuth 
angles <y n &y NUMBER OF 
ARRIVAL DIRECTIONS PER ORBITAL 
STEM (Madpos), AND GEOMAGNETIC 
ACTIVITY LEVEL (Kp AND Dst). 



7 

300 



305 



CHOOSE NUMBER OF STEPS PER 
EACH ORBIT (N spo ) OF THE EARTH 
(E.G. N spo = 200 STEPS), SET ORBIT 
STEP COUNTER N = 0 AND NUMBER 
TRANSMITTED NtrO, AND SPECIFY 
REQUIRED ACCURACY OF NUMERICAL 
INTEGRATION (TYPICALLY 10-5). 



I 



INPUT ORBITAL PARAMETERS: 
INCLINATION.APOGEE, PERIGEE 
INITIAL TIME (to), ORBITAL DURATION (T), 
INITIAL DISPLACEMENT OF PERIGEE FROM 
ASCENDING MODE, INITIAL LONGITUDE OF 
ASCENDING NODE, AND INITIAL DISPLACEMENT 
FROM THE ASCENDING MODE. 



"7 

310 



CALCULATE ORBITAL PERIOD (P) 
FROM INPUT ORBITAL PARAMETERS 
USING STANDARD ORBIT GENERATION 
ROUTINES (E.G. ADAMS ETAL, 
(1986). CALCULATE: 

Nmax = N S po Madpos T/p 

and St = P/N spo . 



315 

V 



320 



FIND SPACECRAFT'S ORBITAL LOCATION 
FOR TIME t=to + N5t USING NUMERICAL 
SOLUTIONS OF STANDARD EQUATIONS FOR 
ORBIT GENERATOR, E.G. THOSE GIVEN ON 
PP. 43-48 OF (STERN, 1960), AS IMPLEMENTED IN THE 
CREME ORBIT GENERATOR (ADAMS ETAL., 
1986). SET PARTICLE TRAJECTORY TIME 
t-t AND TIME STEP 6 1'= 0. SETM = 0. 



DIFFERENTIAL EQUATION FOR INTEGRATING IS 
THE LORENZ EQUATION IN MAGNETIC FIELD, 
USING "BACKWARD" SOLUTION (Q REPLACED BY -Q, 
f REPLACED BY -t') OF: 

F= m 7dv/dt'= Q v x B 
B=B|GRF(r,n + BTsyg89(Kp.r i n 
+ S B EXTENDED Tsyg89(Kp, Dst, r, t') 

NOTE: IN THE "BACKWARD" SOLUTION, 
THE ACTUAL PARTICLE'S FINAL POSITION IS 
THE INITIAL POSITION ALONG THE NUMERICAL 
INTEGRATION OF THE "BACKWARD" PARTICLE 
TRAJECTORY, AND THE "BACKWARD" 
PARTICLE'S INITIAL VELOCITY IS IN THE 
OPPOSITE DIRECTION AS THE ACTUAL 
PARTICLE'S FINAL VELOCITY. 

SEE TEXT FOR DESCRIPTION OF B-FIELD 
CALCULATIONS. 



T 

335 



330 

L 



SET r = PARTICLE'S FINAL POSITION AND v 
= PARTICLE'S FINAL VELOCITY VECTOR. 
THESE ARE DETERMINED FROM THE 
SPACECRAFT'S LATITUDE, LONGITUDE, 
AND ALTITUDE, COMBINED WITH THE 
PARTICLE'S RIGIDITY, FINAL 9, AND 
FINAL (j>. 



•325 



GENERATE PARTICLE'S FINAL 9 AND 
<|> FOR EACH ORBITAL LOCATION, USING 
RANDOM NUMBER TECHNIQUES, I.E. 
x 1 AND X2 ARE TWO DISTINCT OUTPUTS 
FROMARANDOM NUMBER GENERATOR. 

cose = cose max + 

M (cose min - cose max ) 

<t> = <t>min + x 2 (4»max - kin) 
M = M + 1 

NOTE: MANY MODERN COMPUTERS 
HAVE BUILT IN RANDOM NUMBER 
GENERATOR FUNCTIONS. FOR COMPUTERS 
WHICH DO NOT, ONE CAN USE THE SUBROUTINES 
FOR RANDOM NUMBER GENERATION FROM 
(PRESS ETAL, 1992). 



FIG. 3A 
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EVALUATE B-FIELD AND d v/dt* AT LOCATION 
r AND TIME f ALONG THE PARTICLE'S 
TRAJECTORY, I.E. t' IS THE TIME ALONG THE 
PARTICLE'S TRAJECTORY. 

I 

SPECIFICALLY, PERFORM BULIRSCH-STOER 
NUMERICAL INTEGRATION, USING FORMULAE 
& ALGORITHMS AS PRESENTED ON PP. 718 - 725 
OF (PRESS ETAL, 1992). THE OUTPUTS ARE 
THE PARTICLE'S RADIAL POSITION r, VELOCITY 
v, AND TIME t' = t +5 f AFTER NUMERICAL 
INTEGRATION STEP; AND AN ESTIMATE OF 8 f 
FOR THE NEXT PARTICLE-TRAJECTORY 
INTEGRATION TIME-STEP. 

NOTE: 5 1' IS RE-CALCULATED WITHIN EACH 
EXECUTION OF THE ADAPTIVE BULIRSCH-STOER 
INTEGRATION. 



-340 



■345 



EVALUATE IF PARTICLE TRAJECTORY 
ENCOUNTERED SOLID EARTH (FORBIDDEN), 
ENTERED ATMOSPHERE (FORBIDDEN), CROSSED 
MAGNETOSPHERIC BOUNDARY (TRANSMITTED), 
OR STILL WITHIN MAGNETOSPHERE. 



I 
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IF TRANSMITTED, Ni r = N tr +1 
IF FORBIDDEN, Nfb=Nfb+1 
I 



-355 



IF WITHIN MAGNETOSPHERE. 

IF TRANSMITTED OR FORBIDDEN, 
N=N+1 
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IF M < M a dpos 
I 



IF N < N max 
IF N > N max 



I 



-365 



-370 



CALCULATE GEOMAGNETIC 
TRANSMISSION (GT) FOR THIS 
RIGIDITY (R) AND OUTPUT THE 
RESULT: 

GT(R) = Ntr/Nmax 
NOTE:Nfb = N max -Ntr 

FIG. 3B 
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800- 



INPUT: ATOMIC NUMBER (IZ) 




KINETIC ENERGY (En) 




BASELINE MODEL = 


WORST DAY, 




WORST WEEK, 


OR 


PEAK FLUX, 



805 



810 




JZ>2d 



SELECT Fe SPECTRUM 
AS MODEL 



815 



SELECT 0 SPECTRUM 
AS MODEL 



LOOK UP ELEMENTAL 
BREAKPOINT EB2 
FROM TABLE 



-820 



840 

A. 



£ 



845 



LOOK UP A3, 
7-: FROM TABLE 





855 



CALCULATE UNSCALED 
FLUX=A3'EN"Ysi 

1 





LOOK UP A2, 
G s j FROM TABLE 


860 
/ 








CALCULATE UNSCALED FLUX 
=A2'EXP(-G s i'EN 1/4 )'EN 1/4 




FIG. 8 
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B 





A 900 
, / 






LOOK UP A sp , G sp 
FROM TABLE 










CALCULATE UNSCALED FLUX 
=A sp -E N (-Gsp) 



920 

A 




LOOK UP FROM TABLE 
SF (IZ) & SF 
(MODEL ELEMENT) 



930 



925 

7_ 



SOLAR HEAVY ION 
FLUX = UNSCALED FLUX 



CALCULATE SOLAR HEAVY ION FLUX 
+ UNSCALED FLUX - SF (IZ) 



SF (MODEL EL.) 



FIG. 9 
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CREME96 MAIN MENU 

SEE BRIEF INSTRUCTION AT THE BOTTOM OF THIS PAGE. SEE HOW TO RUN CREME96 FOR MORE DETAILS 



60 RESET FORM 


ROUTINE 


USER REQUEST 
FILE LIST 


EDIT USER 
REQUEST FILE 


RUN ROUTINE 


TPS 


▼ 


O 


□ 


GTRN 


▼ 


O 


□ 


FLUX 


T 


o 


□ 


TRANS 


▼ 


o 


□ 


LETSPEC 


▼ 


o 


□ 


PUP 


T 


o 


□ 


HUP 


▼ 


o 


□ 


UTILITIES OT 


\TE PLOTS, PERFORM 


DOWNLOADS, ETC.) 


□ 



1 . TO CREATE A USER REQUEST FILE: 

CLICK ON THE EDIT BUTTON BY THE CORRESPONDING ROUTINE, LEAVING THE USER 
REQUEST FILE BLANK. (IF NECESSARY, USE THE BLANK SPACE IN THE PULLDOWN MENU 
TO CLEAR A NAME ALREADY TYPED THERE.) 
CLICK ON THE 'GO" BAR. 

2. TO EDIT AN EXISTING USER REQUEST FILE: 

CLICK ON THE EDIT BUTTON BY THE CORRESPONDING ROUTINE. 
SELECTAUSER REQUEST FILE FROM THE APPROPRIATE PULL-DOWN MENU. 
CLICK ON THE "GO" BAR. 

3. TO RUN ONE OR MORE ROUTINES: 

IF AN EDIT BUTTON HAS BEEN CLICKED "ON", USE THE RESET BAR TO CLEAR IT. 
SELECT USER REQUEST FILE(S) FROM THE APPROPRIATE PULL-DOWN MENU(S). 
CLICK ON THE RUN BUTTON(S) IN THAT LINE(S). 
CLICK ON THE "GO" BAR. 

4. TO ACCESS CREME96 UTILITIES: 

CLICK ON THE RUN BUTTON ON THAT LINE. 
CLICK ON THE "GO" BAR. 

NOTE: YOU CAN CREATE OR EDIT ONLY ONE USER REQUEST FILE ATA TIME. BUT YOU CAN 
SUBMIT SEVERAL FILES FOR RUNNING (SEQUENTIALLY, IN THE ORDER SHOWN IN THE TABLE) 
BY CHECKING MORE THAN ONE "RUN" BUTTON AND THEN CLICKING ON THE "GO" BAR. 

QUESTIONS? COMMENTS? SEND MAIL TO tvlka@cfs2.nrl.naw.mil. 
YOU MAY TERMINATE YOUR CREME96 SESSION AT ANY TIME. 



FIG. 10 



